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ABSTRACT 



Aims. If NGC 6397 contains a large fraction of "second generation" stars (>70% according to recent analysis), the 
helium abundance of its stars might also be affected, show some star-to-star variation, and be larger than the standard 
Big Bang abundance Y~0.24. Can we derive constraints on this issue from the analysis of the main sequence width 
and from its luminosity function? 

Methods. We build up new models for the turnoff masses and the main sequence down to the hydrogen burning minimum 
mass, adopting two versions of an updated equation of state (EOS) including the OPAL EOS. Models consider different 
initial helium and CNO abundances to cover the range of possible variations between the first and second generation 
stars. We compare the models with the observational main sequence. We also make simulations of the theoretical 
luminosity functions, for different choices of the mass function and of the mixture of first and second generation stars, 
and compare them with the observed luminosity function, by means of the Kolmogorov Smirnov -KS- test. 
Results. The study of the widht of the main sequence at different interval of magnitude is consistent with the hypothesis 
that both generations are present in the cluster. If the CNO increase suggested by spectroscopic observation is taken into 
account the small helium spread of the main sequence in NGC 6397 implies a substantial helium uniformity (AY~0.02) 
between first and second generation stars. The possible spread in helium doubles if an even larger increase of CNO is 
considered. The luminosity function is in any case well consistent with the observed data. 

Key words, stars: evolution - stars: main-sequence 



1. Introduction 

Our views about Globular Clusters (GCs) are dramatically 
changing in these latter years, thanks to precise photo- 
metric investigations, that revealed the pres ence of multi- 
ple main sequences or sub giant branches (e.g. lPiotto et all 
l2007t iMilone et all l2008f ). and the increasing amount of 
new spectroscopic data on GC stars. These data, in par- 
ticular the spectra for about 2000 stars in 19 GCs recently 
obtained by the multiob ject spectrograph FLAMES@VLT 
(|Carretta et all l2009al fbT) have shown how the "chemi- 
cal anomalies" among GC stars are indeed ubiquitous 
in all clusters, and concern a large fraction (from 50 to 
of stars. The peculiar chemical abundances take 



the form of anticorrelations between O and Na, Al and 
Mg, and are present in stars of various evolutionary 
phases (both unev olved main sequence and evolved red gi- 
ant branch stars) ( Gratto n et all 120011 : iRamfrez fc Cohenl 



2002; ICarretta et all [2004) , supporting the idea that they 
are not due to an "in situ" deep mixing in the stars, but 
have been imprinted in the gas from which they formed, 
that was polluted by the winds lost by a previous first 
generation (FG) of stars: according to this hypothesis, we 
are now seeing a second generation (SG) of stars mixed 
to the FG. It is not definitely settled what kind of stars 
produced this material, that must have been processed by 
the hot CNO cycle and other proton capture reactions on 



light nuclei. The two most popular candida tes are massive 
asym ptotic giant branch (AGB) stars (e.g. IVentura et all 
|2001[ ). and possibly , for the extreme anomalies, super- 
AGBs dPumo et all 20081). o r mass ive stars, either fastly 
rotating (e.g. iDecressin etafl 120071). or in binar i es un der- 
going non conservative evolution ~ Jde Mink et al.l . l2009D . 
ICarretta et alj (|2009al ) examine the famous Na-0 anticor- 
relation and tentatively divide the stars of each cluster in 
three groups: the "primordial" stars, those having abun- 
dances of O and Na similar to those found in the halo 
stars; the "intermediate" stars, having high Na and some- 
what depleted O; the "extreme" stars, having high Na and 
strongly depleted O. While only a few, very massive, clus- 
ters contain stars with extreme anomalies, all clusters show 
a population with "intermediate" chemistry. The problem 
of GC formation and early evolution is very complex, but 
it is easy to accept that a cluster contains two or mul- 
tiple populations, if it shows both chemical peculiarities, 
as discussed above, and photometric peculiarities: the GC 
NGC 2808 is a prototype of this clas£|, as it shows three 



1 We do not wish to include to Cen among the classic GCs, 
as it also shows large metallicity variations, indicating that its 
evolution is partially similar to that of a massive GC, containing, 
e.g., a population with a very high helium content, but also 
closer to that of a small galaxy, as the supernova ejecta take 
part in the star formation events. 
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populations both in the main sequence (|Piotto et all |2007| ) 
and in its horizontal bra nch (HB) extended morphology 
(|D'Antona fc CaloU |2004) . and it is possible to reproduce 
these three popula tions by assuming t h at the y differ in he- 
lium content (e.g. iD'Antona fe Caloil . [2008); in addition, 
this cluster shows the Na-0 anticorrelation in one of its 
most extre me forms, with stars r eaching very low oxygen 
abundances. ICarretta et all (|2006f ) even find a possible indi- 
rect hint of helium enhancement in the subgroup of oxygen- 
poor red giants that they examine. As theory expects that 
some helium enrichment accompanies the hot-CNO nucle- 
osynthesis in both the proposed models (massive AGBs or 
massive stars), NGC 2808 represents in many respects, the 
most classic example of well understood multiple popula- 
tions in GCs. 

There are other cases, however, in which the presence 
of subpopulations is not as well clear, and, in particu- 
lar, it is not clear whether the abundance anomalies (in 
so me cases less p r ominen t, but always present, according 
to ICarretta et all l2009a| ) are always accompanied by he- 
lium enrichment, and how large this enrichment is. In par- 
ticular, the cluster NGC 6397 has always been considered 
the perfect example of a "simple stellar population" (SSP) 
due to the "tightness" of its HR diagram, including a very 
compact blue HB. In recent years, spectacular data for the 
low main sequence of NGC 6397 have become available. By 
using the technique of proper motion cleaning. iRicher et alj 
( 2006) obtained a very tight main sequence and a very clean 
luminosity function d own to the hydroge n burning min- 
imum mass (HBMM) (|Richer et all 120081 ). We decided to 
use these data to quantify at what level we can accomodate 
helium variations in this cluster (and the possible, associ- 
ated, CNO variations), by investigating whether the main 
sequence width and its luminosity function are compatible 
with a helium spread and/or with a non-standard helium 
content. At the same time, we take the occasion of this com- 
parison to compute and test new stellar models for the low 
main sequence. In these models, we employ and compare 
new equations of state (EOS) today available, and we test 
the available color-T c g- transformations. 
The outline of the paper is the following. In Section 2 we 
describe in some detail the spectroscopic results concern- 
ing the cluster, and what we expect concerning the multi- 
ple populations it should hide, and what are plausible he- 
lium and C+N+O variat ions expected . We th en summarize 
the photometric data bv lRicher et a l. (2008). After having 
summarized the literature concerning the low mass main 
sequence models, in Section 3 and 4 we describe our code 
and the results of the computation of solar scaled mod- 
els; in Section 5 we present a-enhanced models computed 
for the comparison with the data for the Globular Cluster 
NGC 6397. Both the comparison of the CMD and the lu- 
minosity function with theory are discussed in detail in 
Section 6 also in the hypothesis of multiple populations 
with different helium, and possibly also C+N+O, content. 
In Section 7 we summarize our results and conclusions. 

2. The case of NGC 6397: apparently, a simple 
stellar population 

ICarretta et all (|2009al ) find that at least 70% of stars in 
the cluster NGC 6397 are "intermediate" according to 
their definition, although photometric studies show that 
all the evolutionary sequences in this cluster look like 



those of a prototype SSP: the main s equence is very tight 
(|King et all 119981 : IRicher et all . I2006D . and the horizontal 
branch (HB) lacks the extreme HB and blue hook stars 
that are now regarded as the pro of of the presence of a 
very helium enrich e d pop ulation (jD'Antona et all [2002; 



D'Antona fc Caloil. 2004D . The chemical anomalies of 
Carretta et alj ( 2009al ) analysis, however, do not come as 
a full surprise, as already many hints were available in 
the recent literature about t he dubious simp li city o f this 
cluster. First of all, already [Bonifacio et all ([2002) had 
noticed the presence of Nitrogen rich stars that have a 
nor mal Lithium content (see also Pasquin i et all [2008), 
and ICarretta et all (|2005[ ) find that only three subgiants 
out of 14 stars are nitrogen normal. These features lead us 
to suspect that the material from which these stars formed 
is CNO processed, as expected in the stars with low oxygen 
and high sodium. The possible helium enhancement is 
certainly not extreme, as, for instance, are small li thium 
variation among the turnoff stars (Lind et all 120091 : 
iKorn et all . 120071 : iPasauini et all 120081 ). In the massive 
AGB model for the formation of the second generation, 
the lithium content of the AGB ejecta is not extremely 
depleted as in the other models, but it is difficult to believe 
in a cosmic conspiracy producing exactly the same lithium 
of the FG, unless the AGB matter is very diluted with 
FG gas, so that both lithium and helium do not differ too 
much in the two generations. 

A different, m o stly theoretical, approach, led 
ICaloi fc D'Antona! (|2005h and then ID'Antona fc Caloil 
(|2008l ) to provocatively propose that all the stars in the 
clusters having entirely blue HBs are composed by SG 
stars. This idea is at the basis of a possible explanation 
of the peculiar difference between the GCs M 3 and 
M 13. The cluster M 3 has a complex HB, including 
many stars redder than the RR Lyrae (a red clump), 
RR Lyrae stars, and a well populated blue side, while 
M 13, having the same metallicity, has an only-blue HB. 
This differenc e, that has been generally attribu t ed to 
different age (jJohnson fc Boltd . 11998k iRev et all l200lh 
or to different mass loss along th e red giant branch 
(|Lee et all 11994 ICatelan et all Il998l) (the famous second 
parameter problem) may also be interpreted by assuming 
that M 13 is totally deprived of its FG, and the SG has a 
minimum helium abundance in mass fraction Y~0.28. In 
M 13, whose HB shows a prominent blue tail, simulations 
of the HB stellar distribution show that there must 
also be a small fraction of stars with helium Y~0.38 
(|D'Antona fc Caloil . l2008fh Is it possible that a GC is 
composed only of SG stars? This could indeed happen, as 
shown in some hydrodynamic plus N-body simulations 
of the SG formation and of the cluster first phases o f 
dynamical evolution presented by iD'Ercole et all ([2008). 
They find that the dynamical evolution of the cluster 
may be characterized by expansion of the FG star system 
due to the SNII mass loss preferentially occurring in the 
cluster central regions, while the SG is still forming in the 
core. Depending on the initial conditions, in some cases, 
the ratio of SG to FG stars may even reach a factor 6 
or more. While this kind of modelling depends on the 
input parameters, and does not imply that this really 
occurred in nature, NGC 6397, having a short blue HB, 
and a tight main sequence (MS) and red giant branch 
(RGB), could indeed be made by a homogeneous set of SG 
stars, corresponding to a unique value of Y, just a bitter 
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larger than the Big Bang abundance (Y^O.26-0.28). This 
speculation would help to understand why the nitrogen 
abundance in most of the NGC 6397 stars is quite large. 
The possibility that all the stars in this cluster have a 
homogeneous, but larger than standard, helium abundance 
could be falsified by looking at the main sequence width, 
that depends on the combination of the photometric errors 
with the possible star to star differences in helium. 
Apart fr om this extreme and pr ovocative suggestion, very 
recently ICarretta et al.l (|2009af ) show that only up to 
~30% of stars should belong to the primordial population 
(the FG), in the assumption that all the stars within 
three sigma from the lowest sodium abundance measured 
in a cluster are "primordial'^. So, at present, the most 
reasonable assumption is that NGC 6397 has at least 70% 
of SG. 

Given the sodium and oxygen abundances of the anoma- 
lous cluster stars, do we expect that they have an 
enhanced helium content? In the hypothesis that the SG 
is born from matter mixed with the hot-CNO processed 
e.jecta of massive AGBs, w e can look at the results by 
IVentura fe D'Antonal (|2009f ). They interpret the anoma- 
lous Na-0 abundances in NGC 6397 as a result of mixing 
between 50% of pristine gas with 50% of gas ejected by the 
5M Q AGBs. Looking at Table 2 of their paper, the helium 
abundance in the 5M ejecta for Z=0.0006 is Y=0.329. A 
dilution by 50% with matter having primordial Y=0.24 
provides indeed Y=0.285. If we take these results at face 
value, also the total CNO content of the SG stars is 
larger than that of the FG stars. The 5M Q AGB evolution 
provides a CNO enhancement by a factor ~3, so we must 
consider also increased CNO (and total metallicity) by a 
factor 1.5, when we compute the higher helium models. 
Of course, the computed AGB models do not give a 
mandatory prescription of what really happens in the 
cluster, so that we will also consider normal CNO models 
and models with even larger CNO in order to include all 
possible cases. Helium abundances equal or larger than 
the standard one will be considered, up to Y=0.28, to 
understand whether the hypothesis that at least ~ 70% of 
stars in NGC 6397 have an helium abundance larger than 
the Big Bang abundance is consistent with the photometric 
data. 



2.1. The observational color-magnitude diagram of Globular 
Cluster NGC 6397 

The first HST observations of the low MS of NGC 639 7 
date back to Paresce. De Marchi fe Romaniellol (|1995|) . 



Afterwards, King et alj (| 19981) have observed NGC 6397 



with WFPC2@HST finding that the luminosity func- 
tion has a rapid decline at low mas s end. Recen tly, the 
deepest observations by iRicher et al.l (|2006l . l2008f) refined 
the data. They observed an outer region of NGC 6397 
with ACS@HST using the photometric filters F814W and 
F606W. The high sensitivity of the camera, the large num- 
ber of orbits obtained (=126) and the vicinity of the cluster 



2 Indeed, while ICarretta et al.1 l|2009al ) have only 4 "primor- 
dial" stars in their sample of O-Na measurements for this clus- 
ter, they have determined Na abundances in a larger number of 
stars, and also in this larger sample sodium is "normal", that is 
similar to the sodium of the halo stars having similar metallicity, 
in ~ 25-30% of stars. 



(it is the second closest globular cluster, after M4) made 
it possible to reach the deepest intrinsic luminosities for 
a globular cluster achieved until today, and what appears 
as the termination of the MS. The field observed over- 
laps that of archival WFPC2 data from 1994 and 1997, 
which were used for the proper-motion-cleaning of the data. 
This technique, applied to the deep ACS photometry, pro- 
duces a very narrow MS till its end. These observations 
are a good basis to test the physics of low m ass stars and 
the p ossible role of the helium abundance. IRicher et al.l 
(2008) analyzed the color-magnitude diagram (CMD) dia- 
gram using the results of models comp uted with Dar t mout h 
Stellar Evolution Program (DSEP) (|Dotter et all 120071 ). 
Their main results are that: 

— the MS appears to terminate close to the CMP loca- 
tion of the HBMM predicted by models. The authors 
state that they would have found fainter MS stars in 
the cluster, if there had been any; 

— the MS fitting technique provides a good agreement 
down to F814W=22.5 mag; below this, down to ^24 
mag, the isochrone is either too blue or too low in lu- 
minosity. The authors underline that this is likely due 
to low mass models being less luminous at a given mass 
than real stars; 

— exploring the MS luminosity function, they find that a 
power law for the mass function MF well reproduces 
the distribution in luminosity of the observed stars, 
while a more top-heavy MF is necessary to fit the 
data in the cluster core that they have also available. 
However, theory predicts more stars than observed at 
the lowest MS luminosities, as also p reviou sly found by 
iMontalban. P'Antona fc Mazzitellil (|2000D using data 
bv lKing et al.l (|l998l) ~ 



3. The low mass main sequence models: Input 
Physics of the models 

The computation of very low mass stellar models re- 
quires an accurate knowledge of the EOS for partially 
ionized gas at densities where the ideal gas EOS ap- 
proximations breaks down dramati cally, especially close 
to the pressure ionization region ( Fontai ne et aTT. 119771 : 
Maeni fc Mazzitellil . [l979t IS aumon. Chabrier fc van Hornl 



19951 ). The general properties of very low masses have 
been des cribed in many works (for example see the re- 
views bvlChabrier fc Baraffd . Il997t lAlexander et all 119971: 



ICassisi et al. . 20001 ); for the p o pulatio n II low masses, 



after the work by iP'Antonal (|1987l ). based on grey 
atmosphere boundary co n ditions and o n the EOS by 
iMagni fc Mazzitellil (11979D. iBaraffe et al.l (I1997D presented 
models based on the lSaumon. Chabrier fc van Hornl {l995) 
EOS and on the Ne xtGen non-grey atmosphere model s 
later on published by lHauschildt. Allard fc Baronl ([l999). 
They showed that this latter improvement was essen- 
tial to reproduce the colors of the low mass m ain se- 
quence. Montalb an. P'Antona fc Mazzitellil (|2000f ). reex- 
amining the problem of the EOS, cautioned about the 
use of the Additive Volume interpolation needed to ob- 
tain the thermodynamic quantities for intermediate com- 
positions from the pure hydrogen and pure helium ta- 
bles available in the ISaumon. Chabrier fc van Hornl (fl995l) 
EOS. In recent years, two new EO S have become avail- 
able: the FreeEOS by lITwir] ([2004) and the OPAL EOS, 
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namely the EOS provided by the Liverm ore group as 
a by product of the opacity computation (Ro gers et all 
1996). Both EOS however do not cover the pressure ion- 
iz ation region, for which the best appr oach remains that 
bv lSaumon. Chabrier fc van Horn! (I1995D The FreeEOS has 
been recently employed bv lDotter et al.1 (|2007T ) and applied 
to the fit of the main sequence of NGC 6397. The OPAL 
EOS has not yet been used to approach the construction of 
low mass models, so we decided to adopt it in two different 
ways in our new models, as described in Sect. 4. 
We use the ATON code fo r stellar evolution; a detailed de- 
script ion can be found in IVentura. D'Antona fe M azzitclli 
(2007); in the following we remember the main updated 
inputs that are important for the treatment of low mass 
stars. 



3.1. Opacities and nuclear reaction 



The p rogram includes the opacities by iFerguson et al.1 
(|2005f ) for the external region of the star (T<15000 K) 
and the latest version (2005) of OPAL opacities for higher 
temperatures (|Iglesias fe Rogersl 1 19961 ) . For fully convec- 
tive low mass stars (below ~0.35M Q ), the uncertainties on 
radiative opacities have negligible influence on the models, 
while for larger masses an error in opacity by ~20% may 
cause an erro r up to ~1% in the determination of radii 
(|Dotterl . 120071 ). Electron conduction opacities were taken 
fro m the WEB site of Pot ekhin (2006) and correspond to 
the iPotekhin et~ai1 (|1999( ) treatment, corrected following 
the improvement of the treatment of the e -e scattering 
contribution described in lCassisi et al.l (|2007f ). 
Although not necessary in this computation, the nuclear 
network includes 30 chemical elements, all main reactions 
of p-p, CNO, Ne-Na and Mg-Al chains and the a capture 
of all nuclei up 26 Mg. Th e relevant cro s s sect ion are from 
the NACRE compilation (|Angulo et al.Lll999h . 



3.2. Equation of state: OPAL vs 

\Saumon. Chabrier & van Ho rn (1995) 

As remarked previously, non-ideal effects become increas- 
ingly important for masses M ^ 0.8M©. 
ATON uses 18 tables of EOS in the (gas) pressure- 
temperature plane corresponding to three different metal- 
licities, Z=0, 0.02 and 0.04, and six hydrogen mass fractions 
X, ranging from to 1-Z; the thermodynamic quantities, i.e. 
density, adiabatic gradient, specific heat at constant pres- 
sure C p , the C p /C„ ratio, and the three exponents Fx, T 2 
and T3 are obtained via 4 cubic unidimensional splines on 
X, Z, pressure and temperature. 

These 18 tables are built up in three steps. First, the ther- 
modynamic quantitiejsj^ to the for- 
mulation bv lStolzniarin fe Blocker! ((2000) , which is the most 
modern and updated description available for ionized gas, 
including both classic and relativistic degeneracy, coulom- 
bian effects and exchange interaction. The tables are then 
partially overwritten by the OPAL EOS i n the whole do- 
main where this is available (|Rogers et all Il996l see OPAL 
WEB page, last update in February 2006). Finally in the 
very low-temperature regime, where OPAL EOS is not 
available (see Fig. [H) the tab l es are overwritten by the 
ISaumon. Chabrier fc van Horn! (|1995l ) EOS, which has the 




-5 
Log( p) 



Fig. 1. We show the region of the Log(/?)-Log(T) plane 
(density p in g cm -3 and temperature T in K) where OPAL 
(solid box) and SCH (dashed box) EOS are available. The 
structure profile of two stars with Age=10 Gyr, Z=0.0006 
and with M=0.6 M Q and 0.1 M© are also shown. The 
O.6M0 is in the region of overlapping of the two EOS, while 
0.1M© is exactly at the border between the two EOS. 



advantage of employing an adequate physical model for the 
pressure ionization. The Saumon et al. EOS is only given for 
pure hydrogen and pure helium mixtures; the presence of 
metals is thus simulated by adding helium, and the different 
H-He mixtures must be interpolated through the Additive 
Volume law. 

We call "EOS+OPAL" these tables, that represent 
the standard in our computations. To investigate 
how the results depend on the chosen EOS, we 
built additional tables (EOS+SCH) using the EOS by 
ISaumon. Chabrier fc van Horn! (|l995l SCH) in the whole 
region of the P-p plane for which it is available. This is 
an interesting test, because the structure of the majority 
of stars discussed here are contained in the region of plane 
LogT-Logp where both EOS are available (see Fig. [I]). 
In Fig. [5] we show the differences in the molecular weight (p) 
and adiabatic gradient V a d=dlogP/dlogT along the struc- 
ture of two main sequence models of M=0.3M©, Z=0.0006 
at the age of lOGyr, computed with the two different EOS. 
The differences are more evident in the zone of partially 
ionization where physical differences of the two treatments 
(different p) affect V a d- In the EOS+SCH models, the re- 
gions in which V a d is lower prevail, and the global effect is 
to produce T e g larger by about 100 K, practically at the 
same luminosity, since the inner structure does not change 
significantly with the EOS. 

The results for different masses are shown in Fig. 3, report- 
ing the HR location of the models at lOGyr. The differences 
are small for M>0.5M©, where the regions of partial ion- 
ization do not dominate, and vanish at the lowest masses, 
because only the EOS+SCH is available for the physical 
conditions of their interiors. 

Based on these results, we may conclude that use of both 
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Fig. 2. Adiabatic gradient (V ac i; dash-dotted and dotted 
curves; right y-axis) and molecular weight (/i, solid and 
dashed curves; left ?/-axis) as a function of pressure P (cgs 
units) for the structure with M=0.3 M , Z=0.0006, Y=0.24 
and Age=10 Gyr. The solid and dash dotted curves corre- 
spond to EOS+OPAL while dashed and dotted curves are 
calculated with EOS+SCH (see text). 



g. 3. MR diagram location at 10 Gyr for models with 
Z=0.0006, the labelled masses and calculated with different 
EOS tables. Models from the Dartmouth Stellar Evolution 
Database (full dots) and from Baraffe et al. (1997) (as- 
terisks), with same metallicity of ours, are also shown for 
comparison. 



EOS leads to models with compatible effective tempera- 
tures and luminosities. In FigO for the available masses, 
we report the location of the m odels of same chemistry com- 
puted by iDotter et all (|2007t ) with the Dartmouth Stellar 
Evolution Program (DSEP), using the FreeEOS by Irwin 
(2004) and otherwise a very similar input physics (L and 
T e fi's are taken fro m their WEB p a ge) ; fo r the same masses 
also the models by iBaraffe et al.l (|l997t ) are shown. Note 
that these latter models do not differ significantly from 
ours, both in L and T ff, while the models calculated with 
DSEP at the lowest masses differ, especially in T ff, up to 
^300 K at 0.15M Q . This effect can be possibly attributed 
at the different EOS, but a detailed comparison of models 
would be required. 

3.3. Convection 

The ATON code allows us to model turbulent convec- 
tion by adopting the traditional Mixing Lenght Theory 
(MLT, Bohm-Vitense, 1958) or the "Full Spect r um of 
Turbulence" (FST) model dCanuto fc MazzitelUl . Il99lt 
ICanuto. Goldman &: Mazzitellil . Il996f ) which takes into ac- 
count the full eddies energy distribution (see Canuto e 
Mazzitelli, 1991 for a detailed description of the physical 
differences between the two models). 

While for very low mass stars the description of convection 
has no influence on the atmospheric structure, this does not 
hold for those masses in which convection has a substantial 
degree of overadiabaticity, especially for the stars at the 
turnoff of GCs. For these models a homogeneous modelling 
of convection in the atmosphere and the interior is highly 



recomme nded. At T e ff>400 0K, the grids of models com- 
puted bv lHeiter et~aT1 (|2002l ) by means of NEMO, a modi- 
fied version of Kurucz's code, are available. These grids are 
provided both with MLT model and w ith the FST model by 
ICanuto. Goldman fc M azzitelli (|l996l h A preliminary ver- 
sion o f these latter grids has been used by IMontalban et"aTI 
(|200lD and will be used in this paper (where possible). 

3.4. Atmospheric structure and boundary condition 

At T c ff ;$ 4000K, in the outermost layers of very low 
mass stars, radiative absorption is dominated by molecules, 
and the outcoming flux is far different from the 
frequency-aver a ged d i stribution provided by grey models 
(IBaraffe et all . Il997t iMontalban. D'Antona fc MazziteTiH , 
2000). At larger T e ff, the atmospheric models become less 
critical, whereas the T e ff itself is heavily influenced by the 
treatment of overadiabatic convection. For the non-grey 
models that employ an MLT treatment of convection in the 
atmosphere with a given a a tm , not only the a = / / H p used 
in the interior computation (ai n ) affects the T e ff, but also 
the optical depth at which the match between the atmo- 
spheric a nd the interior in t egrati o n is made, and th e value 
of a atm (|Montalban et all I2004D . IMontalban et alj (|200lh 
have sh own that use of the NEMO grids of model atmo- 
spheres (|Heiter et all I2002T ) computed with the FST con- 
vection may provide a good match to the interior models 
computed with the same convection model in the interior, 
independently of the matching optical depth. 

For the above two reasons we use, according to the stel- 
lar mass, boundary conditions based on two different grids 
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of non grey models of atmospher^B For M> 0.5M Q we use 
iHeiter et all (|2002t ) FST grids, and FST convection also 
in the interior computations, whereas for M< 0.5M(7) we 
adopt the NextGen grids by [Hauschil dt. Allard fe Baro"nl 
(|l999l ) computed by the PHOENIX code with the MLT 
treatment of convection and a a tm — I /Hp = 1. The avail- 
able grids extend down to T off =800 K for the [M/H]=- 
2.0 models, but just down to 2000 K for larger metal- 
licity. For these M< 0.5M© models, we adopt MLT con- 
vection also in the interior computation, setting a in =2.cQ. 
Since the model atmospheres are computed assuming an 
ideal gas and a Saha-like thermodynamics, they should 
not be used in the real gas domain where pressure effects 
are relevant, and a relatively small value of the optical 
depth is chosen for the match between atmosphere and inte- 
rior. We use Tph =3(l0) for M<(>) O.5M . Notice that the 
grids of model atmospheres available are computed only 
for solar scaled mixtures. The lack of suitable model atmo- 
spheres for a-enhanced populations forces us to use for the 
boundary atmospheric conditions a grid with larger [Fe/H], 
to simulat e the a-enhancemen t, following the procedure 
adopted bv lBaraffe et"aLl (j 19971 ). We use a grid for [Fe/H]=- 
1.7, obtained through interpolation between the grids for 
[Fe/H]=-2.0 and -1.5, to compute a-enhanced models with 
Z=0.0002. 



3.5. Transformations to observational plane 

In order to compare the models to the photometric data in 
NGC 6397 we convert luminosity, T e ff and surface gravity 
into absolute magnitudes and colors in the ACS filters. The 
method most used is to calculate theoretical stellar spectra 
from atmosphere's models and convolving these synthetic 
spectra with the filter transmission curves for a photomet- 
ric system which defines the transmission of light through 
the filter as a function of wavelength. Uncertainties de- 
rive especially from the missing or incorrect absorption fea- 
tures and simplifying physical laws, such the assumption of 
LTE. On the other hand, semiempirical colors a nd bolomet- 
ric corrections (e.g. IVandenberg fc Cleml . l2003f ) have other 
uncertainties, e.g. t hey depend on the assumed distance 
potter et all 120071) . 

Fo r the specific case of ACS filters, we use the procedure 
by iBedin et al.l (|2005[ ). They computed ACS color indices 
by using homogeneous set of ODFNEW model atmospheres 
and synthetic fluxe s comp uted with Kurucz ATLAS9 code 
(|Castelli fc Kuruczl . 120041) . Grids of magnitudes for differ- 
ent values of [Fe/H], for 3500K < T cff < 50000K, < logg < 
5.0 and microturbolent velocity 2.0 km/s -1 are provided. 
Visual bolometric correction BCy, visual magnitude My, 
and color indices Mv-Macs ar e given. The ACS magni- 
tudes wer e computed by using the WFC/ACS transmission 
curves bv lSirianni et all ([20 05') . while they adopted the V 
passband from iBessell (|1990l) . Finally, they assumed that 
Vega ACS magnitude would be equal to 0.00 in all pass- 
bands. 



The same approach has been adopted for Pre M ain Sequence 
stars in lDi Criscienzo. Ventura fc D'An tona (20081) 

4 The choice of the a parameter is less and less critical when 
decreasing the mass, as the external layers become so dense that 
convection becomes more and more adiabatic. Nevertheless, we 
use ai n =2.0 to allow for a smooth T e ft transition between the 
upper (M> O.5M ) and lower (M< O.5M ) MS models. 
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Fig. 4. Evolutionary tracks, starting from the pre-main- 
sequence, for three different masses with Z=0.0006 and dif- 
ferent value of initial helium. The filled circle label the po- 
sition in HR diagram of the star at Age=10Gyr. Only for 
the highest mass the tracks with the all three Y are shown 



The bolometric corrections bv lBedin et al. (2005) for ACS 
filters do not extend below T e g=3500K. For these low tem- 
peratures we can use the values obtained from the synthetic 
spectra of lHauschildt. Allard fc Baron (I1999D. Also i n this 
case, the transmission filters of Sirianni et all (120051) were 
used, and zero points corrections to standard system are 
obtained from observed Vega spectra. 

Obviously the availability of a unique set of colors- Teff 
transformations would be highly recommended; however, 
at least, for masses below 0.5 M Q , we have the advantage 
of using the same bolometric corrections derived from the 
atmospheric structures used as boundary conditions for the 
stellar models. 

Finally we note that the two sets of correlations match very 
well in the main sequence, so that no discontinuity arises 
in our transformed isochrones. 



4. Results: solar scaled models 

We computed evolutionary tracks for low mass objects 
(M<O.8M ) from the pre main sequence to the red giant 
branch, or until they reach an age of 20Gyr. Results are 
shown here for [M/H]=-1.5 and [M/H]=-2.0, but larger 
metallicities ([M/H]=-1.00 and -0.50) are available upon 
request. 

In thi s Section we present re sults for solar scaled mix- 
tures (jGrevesse fc Sauvall . [T99l GS1999), [a/Fe]=0.4 mod- 
els will be used in the Section 6 to compare with the data 
of NGC 6397. Models are extended down to the HBMM 
when the atmospheric boundary conditions allow it. 
We computed standard evolutionary tracks with an initial 
helium c ontent close to t he Big Bang abundance (that is 
Y=0.24. ICoc et all l2004f) . and models with larger helium 
(Y=0.28, 0.32 and 0.40). From the evolutionary tracks, 
isochrones are derived for typical ages expected for GCs, 



Di Criscienzo et al.: The helium content in the Globular Cluster of NGC 6397 



7 



_ 1 1 1 1 1 1 1 1 1 1 
L[M/H]=-1.5 

. >* 


i | , , , | , , , _ 
■ C », 


— 1 1 1 1 1 1 1 h 

L[M/H] = -2.0 


• — | — | — | — | — — | — | — 










- 













-1.2 -1 -0.8 -0.6 -0.4 -0.2 

Log(M/M Q ) 

Fig. 5. The MLR for models at 10 Gyr having the la- 
belled value of metallicity and different helium abundance 
(=0.24-open circles, =0.28-filled triangles, =0.32-filled cir- 
cles, =0.40-asterisks). In each case the lower masses are the 
derived HBMM and reported in Table 1. 
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Fig. 6. Comparison at very low masses between the MLR 
obtained using models computed with different equations 
of state (see text). 



from 10 to 14 Gyr. 

Fig. [U compares the HR diagram evolution for three differ- 
ent masses (M=0.70, 0.30 and 0.10 M Q ) and different Y. 
Models with a larger helium abundance have larger lumi- 
nosity and T e ff, due to the average larger mean molecular 
weight [i. This effect is more evident in the stars with a 
radiative core (here shown is the M=0.7M©) than in the 
totally convective stars like the O.3M . The difference in- 
creases again in the lowest masses ( ,$ 0.15 M©), where par- 
tial degeneracy begins playing a role. 

4.1. Mass-luminosity relation 

The mass-luminosity relation (MLR) is essential for the 
comparison with the data, as it enters in the conversion 
of the (assumed) mass function into the luminosity func- 
tion (LF). 

Any change of slope of the MLR will be reflected in the 
luminosity function, defined as dN/dL=dN/dM x dM/dL. 
Where the MLR presents an inflection point, the LF has a 
relative maximum or minimum. As the luminosity decreases 
along the MS, there are two inflection points, responsible for 
two main peaks in the LF (see Figj7j) : the first, at Mpgi4w 
~ 6 mag, is due to the transition between pure MS mod- 
els and models that suffer the effects of evolution (and are 
more luminous in the MS due to the hydrogen consump- 
tion). Therefore, the corresponding peak in the LF is a 
function of age, and, we will see, of Y. After a small range 
of homologous models, the MLR steepens progressively due 
to the onset of molecular absorption in the stellar envelope; 
this produces a gradual increase in the LF. When models 
become fully convective, at ~0.35M Q , the MLR relation 
begins to flatten again, and the presence of this inflection 
point results in the large peak at Mp si4w ~ 8 mag, present 
in all the GC LFs (jD'Antonal . Il998f ) and shown in Fig. 
Finally, at masses M ~0.12M Q , the MLR relation flattens 
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Fig. 7. Theoretically predicted luminosity function for the 
models at Agc=13 Gyr and with Z=0.0002 and two differ- 
ent helium abundance. A power law with exponent a=-0.5 
was used as a mass function. Also the LF at Age=10 and 
with Y=0.24 is shown for comparison. 



even more. This is due to the on set of degeneracy in the 
core of the star (jD'Antonal . 119981 ). and this latter decrease 
in the LF is dependent on the EOS. The HR diagram does 
not show the minute features of the MLR derivative, but 
it shows two "kinks", the first one corresponding to the 
onset of molecular hyd rogen dissociation in the envelope 
(|Copeland et all Il970t ) and the second to the onset of de- 
generacy. 
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Table 1. Hydrogen burning minimum mass (HBMM) for 
Z=0.0002 and Z=0.0006, [a/Fe]=0, and different Y. 





Y=0.24 


Y=0.28 


Y=0.32 


Y=0.40 


Z=0.0002 
Z=0.0006 


O.O88M 
0.085M© 


0.083M© 
O.O81M 


0.079M Q 
O.O77M 


0.072M Q 
0.069M Q 



In Fig. [5] we show the dependence of the MLR on the metal- 
licity and helium content at 10 Gyr. Y influences mostly 
the evolved part of MS and the location of the very low 
masses, where degeneracy sets in and the second MS kink 
is located. At the MS end, decreasing the mass, the higher 
helium models remain more luminous and hotter. These 
features will affect the low luminosity LF, that will decrease 
more slowly with decreasing luminosity for larger Y. At the 
largest MS luminosities, the larger is Y, the smaller is the 
slope of the MLR of those models that partially burn their 
hydrogen during a Hubble time, with respect to the case 
of Y=0.24. This is an obvious feature of the evolutionary 
models: the larger Y produces larger MS luminosities and 
faster MS evolution. 

The quantitative results concerning the lowest luminosi- 
ties will also depend on the EOS, as we can easily un- 
derstand by comparing the MLR relations obtained using 
EOS+OPAL and EOS+SCH tables (Fig. For masses 
larger than O.IMq the models using EOS+SCH are slightly 
more luminous, as a consequence of the differences in V a( j, 
but the trend is reversed closer to degeneracy. Obviously 
this variation in the slope of the MLR will produce a dif- 
ferent shape in the peak of the LF, and in particular we 
expect a stronger peak when EOS+SCH are used. 
An interesting theoretical feature of the models at the 
boundary between low mass stars and brown dwarfs is 
shown in Fig. [51 This figure refers to the Z=0.0002 case, 
for which an extended atmospheric grid is available. In be- 
tween the HBMM — the smallest mass that stabilizes in 
MS, in a configuration in which the total luminosity is pro- 
vided by the proton-proton (p-p) reactions in the core — 
and the pure brown dwarfs — that never ignite the p-p 
chain — there is a small range of masses for which nuclear 
burning contributes to the stellar luminosity for even sev- 
eral billion years, but in the end these objects finally cool as 
brown dwarfs. This is a common o ccurrence in population 
I (the transition masses defined in lDAntona fc M azzitelli. 
fl985l) . where the MS merges without discontinuities into 
the brown dwarf cooling sequences. On the contrary, the 
MS of the population II has a much sharper drop, because 
the much smaller opacities put the HBMM at minimum lu- 
minosity a factor ^10 larger, so that the transition masses 
cover a very small mass range. There is then a "luminos- 
ity gap" between the end of the MS and the luminosity 
at which the smaller brown dwarfs are able to slow their 
cooling down to the typical age of population II stars (10- 
12 Gyr). Fig. [6j translates into possible LF. In Fig[7] LFs 
for 10 and 13 Gyr are plotted, assuming a power law MF 
with exponent a=-0.5. At 13 Gyr, the low mass brown 
dwarfs should emerge as a small peak at Mp8i4W ~18, cor- 
responding to near infrared magnitudes ~30 at the distance 
of NGC 6397. The peak is ~1 mag brighter if the age is 3 
Gyr smaller. We regard this prediction as an educated guess 
on the possibility that the dimmest luminosities regime is 
populated not only by white dwarfs (Richer et al. 2008) 
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Fig. 8. We show the comparison between the diagram to 
the standard CNO tracks at different masses(full line) and 
the "peculiar" ones calculated (dashed and dotted line) 
with the evolution of the same mass. 

but also by very cool brown dwarfs. However, dynamical 
models indicate that these objects should be preferentially 
stripped from the cluster. Since there is a strong difference 
in these two populations, as the cool white dwarfs would 
be located at a color Mf606W-Mfsi4w ~l-1.2mag, while 
the cool brown dwarfs will not be visible in the F606W 
band, having Mf606W-Mfsi4w ~5.5-7mag only future ob- 
servations with even more capable telescopes should clarify 
the question. 

5. a-enhanced models, FG and SG populations, 
CNO enrichment 

Comparison with low metallicity cluster stars requires use 
of models with a-enhanced mixtures, for which we adopt 
[a/Fe]=0.4. As discussed only up to ^30% of stars of 
NGC 6397 can belong to the FG, the majority show in- 
deed the Na-0 anticorrelation, and many stars have very 
high Nitrogen content. Therefore we need also models to 
represent the SG. We will assume that it may differ both in 
helium content and in total CNO content from the standard 
FG. A resonable assumption for CNO is that we assume an 
overbundace of N by about 1.4 dex, and a variation of -0.2 
dex for O, leaving carbon unchanged (A. Bragaglia, private 
comunication). The total "metallicity" Z in mass fraction 
for this CNO-enhanced mixture (CNOT) is now Z=0.0003. 
We used the OPAL Web tool to compute on purpose radia- 
tive opacities for this mixture. For T< T5000K we still use 
the opacities by iFerguson et al.l (|2005l ) as lower tempera- 
ture opacities do not affect the structure of the models we 
are considering. We compute on purpose radiative opacities 
for this mixture. All models are computed for helium mass 
fraction Y=0.24 and Y=0.28 and summariz ed in Table 2. 
In gen eral, as discussed in the analysis by IVentura et al. I 
(2009]), the largest differences between standard and CNO- 
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Name 


Description 


Z 


[C/Fe] 


[N/Fc] 


[O/Fej 


Reference 


Z=0.0006 


CNOxls 


Solar- scaled 


0.0006 


0.00 


0.00 


0.00 


GS1999 


CNOxla 


[a/Fe]=0.4 


0.0006 


0.00 


0.00 


0.40 


GS1999 


Z=0.0002 


CNOxls 


Solar- scaled 


0.0002 


0.00 


0.00 


0.00 


GS1999 


CNOxla 


[a/Fe]=0.40 


0.0002 


0.00 


0.00 


0.40 


GS1999 


CNO| 


total(CNO) = 1.6 


0.0003 


0.00 


1.40 


0.20 


this work 




Fig. 9. Left panel Main Sequence of NGC 6397 from iRicher et all ((2008). Also shown are the best fit isochrones for 
FG (red) and SG (green) populations. Numbers indicate stellar masses (up-FG,down-SG) of same models (crosses-FG, 
open circle-SG) in M Q units. Right panel: NGC 6397 MS Luminosity Functio n from data; for three low luminosities are 
reported the completeness fraction taken from Table 4 of lRicher etail (|2008l ). 

6. Comparison with the data of NGC 6397 

We analyse both the CMD diagram and the luminosity 
function simulations. 

6.1. Color-Magnitude diagram 

In Fig. [5] we compar e the models with t he deep photom- 
etry of NGC 6397 bv IRicher et "all ((20081) . We plotted the 
isochrone that allows the best fit of both the low main se- 
quence and the TO for the labelled value of distance modu- 
lus an d reddening, with [Fe/H]=-1.99 dex ((Carretta et all 
2009c) and [a/Fe]=0.4 (solid line-red in the electronic 
version) and with an age of 12 Gyr which is compara- 
ble with the age obtained from the Whit e Dwa rf Cooling 
Sequence of NGC 6397 bv lHansen et~ai1 (|2007h . The dis- 
tance modulus in F814W photometric band is compat- 



enhanced mixtures are found in the ionization zone of the 
CNO elements, but in this case the differences are very 
small since the variation is not very large, and the initial 
abundances are very small. As a result, pratically no differ- 
ences are found in effective temperature and luminosity of 
the models (see FigE}. 

Another possible hypothesis is that Oxygen in the SG is 
basi cally not very d ifferent from the FG value, as most of 
the ICarretta et al.l (|2009a) measurements for O are only 
upper limits. In this case, the CNO abundance becomes 
even larger. In this case, we adopt as boundary conditions 
the grid [M/H]=-1.5. (CNOTT models, see Table 2). As re- 
ported in Fig[8] also in this case the differences are very 
small but for masses of about 0.3M© the temperature are 
little smaller. 
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ible with the true distance modulus (=12.03 mag) and 
E(F606W-F814W)(=0 20 mag) reported in literature (see 
for example Table 3 of Ri cher et al.l, 120081) . 
As already found by iRicher et alJ (|2008t ). below 0.20 M Q 
the isochrone and the data do not match perfectly, how- 
ever the data appear to terminate at about the magnitude 
predicted by mode ls. We then confirm the suggestion of 
IRicher et al.l (|2008f ) that they have observed the termina- 
tion of the hydrogen burning sequence. 
In Fig. [9] we also show (dashed line-green in the electronic 
version) the position of the isochrone obtained with mod- 
els computed with Y=0.28 and CNOj using the same dis- 
tance modulus and age used for the FG isochrone. The two 
isochrones are very similar except around mF8i4w=23 mag 
where they deviate and the FG isochrone is a little brighter. 
This aspect is important to understand if the tightness of 
the MS at this interval of magnitude depends on obser- 
vational error only or is the consequence of the presence 
of a second generation made of stars with higher helium 
abundance and CNO enhancement. We select data within 
three different intervals of half magnitude below the MS 
turnoff, and rectify their colors by subtracting the color of 
their best-fit line, as shown in Fig. [TQ1 left and medium 
panels. We do not consid er the possible presence of bina- 
ries, as lDavis et"ail (|2008ft have shown that NGC 6397 has 
a primordial binary fraction of only The histogram 

of the color displacements from the best fit line, shown in 
the right panels, are fitted with a gaussian profile with the 
labelled a b s . 

In the same panel, for each interval of magnitude we report 
the color displacements for synthetic populations (assum- 
ing the distance modulus for NGC 6397 reported in Fig[5] 
) under the hypothesis that the 30 % of all stars are pri- 
mordial (standard helium and CNO abundance), and the 
restant 70 % is composed by a mixture of SG population 
with CNOt and helium respectively up to Y=YUP=0.25, 
0.26 and 0.28. As expected, the largest difference between 
a is obtained in the third interval of magnitude. From the 
values of dispersion reported we obtain that a SG made of 
stars with an helium dispersion of AY=0.02 is compatible 
with the tightness of the MS of NGC 6397. 
In the case of CNOft models (see Section 5) the isochrone 
calculated for Y=0.28 overlaps exactly on the FG's one. In 
this case the observational spread of MS at low magnitude 
may suggest an even larger spread of helium between FG 
and SG. 



6.2. Luminosity function 

We now compare the observed luminosity function to the 
theoretical simulations. We use the magnitudes in F814W 
filter and derive synthetic populations from models at dif- 
ferent ages, in the plausible range from 9 to 14 Gyr. 
The mass-luminosity relation from our models was dis- 
cussed in Section 4.1; now we discuss the choice of the mass 
function (MF). At the beginning we consider this cluster 
as formed by an unique population in order to define a 
method to compare theoretical and observational luminos- 
ity function. Due to the dynamical evolution of the cluster, 
the present M F is not the initial one. NGC 6397 has a 
collapsed core (jDiorgovski fc Ki ng. 1986), as it has evolved 
past the potentially catastrophic phase of core collapse, and 
is dynamically old. In particular it has been shown that 
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Fig. 10. a) In each row a little portion of the MS of NGC 
6397 is shown on the leftmost side, while in the middle 
panel the magnitude mpsuw is reported versus the color 
distance from the best-fit line. On the right panel we show 
the histogram of the color displacements, fitted with the 
best fit gaussian curve for the data. We also report the 
a of the best fit gaussian curve obtained for of synthetic 
population made of 30 % of stars belonging to FG and the 
restant 70 % made of stars of SG with different YUP. 



NGC 6397 exhib its mass segregatio n , which certainly has 
affected the MF (jHurlev et al.l . I2008D . Isiivestri et all (|1998T> 
ha ve shown, adopting t heir own low mass models and those 
bv lBaraffe et al.1 (|l997t ). that the bulk of the MF can be de- 
scribed by a unique power law of t he fo rm dN/dM=k M a 
with index a=-0.5. IRicher et alJ ((2008) found that a=- 
0.13 gives the highest \ 2 when comparing the models with 
their data of NGC 6397. They found an even better result 
when a lognormal distribution is used, which have the ad- 
vantage to truncate the LF at the extreme low mass end 
but the disadvantage to introduce a new parameter. A sim- 
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Fig. 11. Distribution of KS with distance modulus for the 
labelled value of Z and Y and for different an ((upper- 
panel)), different ages (medium-panel) and MF8i4W.cut 
(down-panel) (see text). In all cases a double power law 
was used with oe,=— 0.1. In the medium and down panels 
o>h=-0-7 as inferred by the upper panel is used. 



i lar re sult was obtained by iDe Marchi. Paresce fc Pulond 
(2000)which suggested that no single power law distribu- 
tion is compatible with the MS of NGC 6397; in particular 
they found that the MF is less steep for M<0.3M Q . 
We simulate the synthetic populations using a power law 
MF with two different slopes above (an) and below 
a cutoff magnitude in the range MF8i4w,cut=8.5-10.5 mag 
(here we consider a single MF as a particular case with 
an =oti,.); this cutoff magnitude, for the chosen distance 
modulus, corresponds to a "cutoff mass". At a fixed age, 
MLR gives the magnitudes of each extracted mass, accord- 
ing to the chosen MF, which depends by four parameter 
(an, ctL, Mp8i4w.cut and distance ). We simulate photomet- 
ric errors considering gaussian errors for the magnitudes of 
each extracted mass in order to reproduce the width of the 
upper main sequence. In addition we have taken into con- 
sideration the uncertainties in completeness, by multiply- 
ing the random extractions in a given interval of masses for 
the co mpleteness fractions given in Table 4 of R icher et al.l 
( 200 81) and determ i ned w ith the artificial star test described 



And erson et~aT1 (|2008). 



We have then used the Kolmogorov-Smirnov test to com- 
pare the observed luminosity function of NGC 6397 with 
the theoretical one, which depends on the MF used to ex- 
tract masses (an ,o;iMF8i4w,cut), on age and distance mod- 
ulus. This statistical method has the advantage of being 
non-parametric and without making assumptions about the 
distribution function of the data, it returns the probability 
that two arrays of data values are drawn from the same 
distribution. The scalar KS3 (varying between and 1) 
which give the significance level of the KS statistic i.e. the 
probability with which we can accept the null hypothesis. 
KS=1 means that the simulated and observed data follow 
the same function. 

Another great advantage of this statistical method is that 
the two arrays of data does not need to have the same 
number of elements; this means that we can build our syn- 
thetic populations with a greater number of stars then the 
real stars from which we made the observed LF, making 
the results independent from the random extraction. We 
have studied the dependence of KS numbers as a function 
of distance modulus for each of the four free parameter 
(age, an, ctL and Mpsi4w.cut) used to build the synthetic 
population and then we have compared the results in or- 
der to choice the best fit parameter. In particular as done 
by Richer et al. (2008) for their x 2 method both an and 
oll were allowed to range from -1 and 1 in 200 steps con- 
sidering. We have explored the case of a single MF as a 
particular case of this situation (ctH=ctL)- We find that the 
best combination of power law exponents for, respectively, 
higher and lower masses, are a#=-0.7 and «l=-0.1. The 
upper panel of Fig. 11, where for example are shown the 
variation of KS with distance modulus for different selected 
Qg, justify our choice. We also note that the distance mod- 
ulus for which we have the higher probability is about the 
same distance modulus for which we have the best compar- 
ison between isochrones and CMD (see FiglH]), confirming 
the validity of our method. The same result is obtained us- 
ing models with higher helium abundance. For these values, 
in Fig. [Tl]we also report KS as a function of distance mod- 
ulus for different ages (medium panel) and for different cut- 



5 We use the kstwo algorithm from "Numerical Recipes" ,3nd 
edition. 
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Fig. 12. Comparison between observed and calculated dis- Fig. 13. Distribution of KS with distance modulus for dis- 
tribution of MS stars of NGC 6397 for the labelled values ferent synthetic populations, 
of the parameters. The power law used is the same of Fig. 

[TT1 The error bars take into account the Poisson's error and . , . , ,,, ,._ , , ,. ..... 

a mixture ol stars with different helium, or a single helium 
incompleteness corrections. „ , ' to 

abundance for all stars. 



off magnitudes (lower panel). In the bottom panel, also the 
case of a single power law with index -0.7 is shown (black 
triangles): we see that a unique power law does not give a 
good match to the observed LF. Notice however that for 
higher m asses this value is ve ry different from the one ob- 
tained bv iRicher etaT] (|2008t) (=-0.13). However we want 
too stress that this our result is consistent with their con- 
sideration that a lognormal function produces a better \ 2 
values than the best fitting single power law MF. In fact 
when a lognormal function is used one has two adjustable 
parameters as in our case. Concerning the age, as shown in 
medium panel of Fig. 11, the best agreement with obser- 
vations is obtained for an age of 14-13 Gyr, but if we also 
consider the best distance moduli and take into account the 
CMD we can choose 12 Gyr together with Mp8i4w,cut=9.5 
mag, corresponding to Af~O.18M(r0- 

In Fig. [H] we report the NGC 6397 MS LF compared with 
the best fitting double power law mass functions. The over- 
all agreement is satisfactory. Only at MF8i4>12mag (~ 
M=0.090M Q ), the observed stars are fewer then predicted; 
this may mean that the MF is even flatter at these lowest 
masses, due to more effective evaporation, but it may also 
be due to some deficiency in the models. In Fig. [T2] we 
also show the distribution of KS in the case that 30% of 
stars belong to FG (Y=0.24 and CNOxla) and 70% to the 
SG (YUP=0.26 and CNOt) which give the best match of 
the widht of the MS (see Fig. [TO]) and in the case that all 
stars belong to SG. The conclusion is that, since compa- 
rable KS, and reliable distance moduli are obtained in all 
cases, as from the widht of MS also from the luminosity 
function we cannot exclude that the cluster contains cither 



6 This value is much lower then th e mass (~O.3M0) found by 
iDe Marchi. Paresce fc Pulond (12000ft 



7. Conclusions 

We have computed new models for the main sequence 
down to the hydrogen burning minimum mass, adopting 
two different version of an updated equation of state and 
made simulations of the luminosity functions for different 
choices of the mass function and the initial helium content. 
The results are compared with the recent observations 
of the MS of NGC 6397 by Richer et al. (2008). Using 
a Kolmogorov Smirnov test to compare observed and 
simulated LF we found that a double power law for the 
mass function well reproduces the observed luminosity 
function in F814W photometric band. However, both 
the models for simple or mixed population according the 
spectroscopic data provide a good fit of LF. A stronger 
results is obtained from the analysis of the width of MS 
from which we find that in any case any helium variations 
must be confined within AY=0.02 in the case of CNO 
overabundance predicted by a mixing between 50% of 
pristine gas and 50% of gas eject by 5M AGB stars as 
suggested by Ventura & D'Antona (2009). Instead we find 
that a arger spread ( 0.02<AY<0.04) in helium between 
primordial and intermediate generation is compatible with 
the widht of main sequence when CNOff models are 
considered. 

The complete sets of isochrones transformed for ACS fil- 
ters F814W and F606W, calculated for this work, are avail- 
able upon request to the authors and will be soon inserted 
on WEB at site |http://www.mporzio.astro.it/% 7Etsa/ 
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